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ABSTRACT 



5-H 



Context. The detection of gamma rays in the very-high-energy (VHE) energy range (100 GeV-100 TeV) provides a direct view of the parent 
C$ ' population of ultra-relativistic particles found in astrophysical sources. For this reason, VHE gamma rays are useful for understanding the 
underlying astrophysical processes in non-thermal sources. 

Aims. We investigate unidentified VHE gamma-ray sources that have been discovered with HESS in the most sensitive blind survey of the 
Galactic plane at VHE energies conducted so far. 

Methods. The HESS array of imaging atmospheric Cherenkov telescopes (IACTs) has a high sensitivity compared with previous instruments 
(~ 0.01 Crab in 25 hours observation time for a 5<x point-source detection), and with its large field of view, is well suited for scan-based 
observations. The on-going HESS survey of the inner Galaxy has revealed a large number of new VHE sources, and for each we attempt to 
associate the VHE emission with multi-wavelength data in the radio through X-ray wavebands. 

Results. For each of the eight unidentified VHE sources considered here, we present the energy spectra and sky maps of the sources and their 
environment. The VHE morphology is compared with available multi-wavelength data (mainly radio and X-rays). No plausible counterparts 
are found. 
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1. Introduction 
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VHE gamma-ray astronomy has recently entered an new era of 
discovery with the introduction of the latest generation Imaging 
Atmospheric Cherenkov Telescopes (IACTs) such as HESS 
(the High Energy Stereoscopic System). Since HESS began 



2 



The HESS Collaboration: Unidentified HESS Sources 



operation in 2004, about two dozen new VHE sources have 
been revealed. Presently identified VHE gamma-ray sources 
belong to one of four categories: active galactic nuclei (AGN), 
pulsar wind nebulae (PWN), shell-type supernova remnants 
(SNR), or X-ray binaries (XRB); recently also VHE emis- 
sion was detected which may be associated with a young 
stellar cluster (Aharonian et al. 2007c). All these identified 
source classes also exhibit emission in the radio and/or X-ray 
regime. However, several VHE sources discovered by HESS 
in the field-of-view of other known sources (Aharonian et al. 
2005c) or during the HESS Galactic plane survey (Aharonian 
et al. 2005b, 2006d) have not been identified with objects 
from which VHE emission is expected. The first unidentified 
VHE source was TeV J2032+4130 (Aharonian et al. 2002, 
2005a), which was discovered by the HEGRA IACT sys- 
tem. HESS J1303-631 (Aharonian et al. 2005c) was found 
in the field-of-view of the binary pulsar system PSR B1259- 
63/SS 2883, and several other sources were subsequently dis- 
covered in the Galactic plane survey. To date, these objects re- 
main unidentified; HESS J 1303-631 has even been postulated 
to be related to such an exotic phenomenon as a gamma-ray 
burst remnant (Atoyan et al. 2006). 

VHE gamma rays are tracers of non-thermal particle accel- 
eration, and their production can be explained by the presence 
of either high-energy electrons or protons. In electron scenar- 
ios, gamma rays are primarily produced by inverse-Compton 
up-scattering of background photon fields by high-energy elec- 
trons. Significant X-ray and radio emission is predicted since 
the same population of electrons should emit synchrotron ra- 
diation at longer wavelengths. For typical Galactic magnetic 
field strengths, the energy flux of the X-ray component of the 
photon spectrum in the keV range is predicted to be compara- 
ble to the energy flux in the TeV range. The X-ray component 
of the spectrum may be suppressed, however, if there is a cut- 
off in the parent electron spectrum below ~10 TeV (Aharonian 
et al. 1997). In proton scenarios, VHE gamma rays are pro- 
duced primarily from the decay of neutral pions (n°) that result 
from proton-proton interactions. If gamma rays are produced 
only via n° decay, a strong X-ray or radio signal may not be 
present; however, proton interactions also produce charged pi- 
ons and cascades of secondary electrons that should generate 
a continuum of X-ray and radio synchrotron emission. Since 
it is difficult to explain VHE gamma-ray emission without at 
least a weak lower-energy counterpart, the lack of low-energy 
emission from the unidentified HESS sources puts significant 
constraints on physical conditions and/or particle acceleration 
processes in their sources. While the explanation may simply 
be that sufficiently deep multi-wavelength observations of the 
objects have not yet been made, the possibility exists that there 
is a new class of object that does not follow the predictions of 
standard emission models. 

Recent observations of the Galactic plane and further re- 
observations of known sources with HESS have allowed for 
the study of some of the weaker Galactic sources at increased 
sensitivity and have revealed new VHE gamma-ray sources in 
addition to those described by Aharonian et al. (2006d). Similar 
to the previously mentioned objects, several of these sources 
have no obvious cataloged counterpart at longer wavelengths, 



and consequently their emission mechanism is unidentified. 
In this paper, we focus on eight VHE emitters without obvi- 
ous counterpart that have been detected by HESS. Of these 
sources, an updated analysis is given for two previously pub- 
lished unidentified sources for which subsequent observations 
have provided significantly better statistics, and the detections 
of six new unidentified sources are reported. New VHE detec- 
tions within the Galactic plane of known objects (PWN, SNRs, 
etc.) have been or will be reported elsewhere (e.g. in Aharonian 
et al. 2007a,b,c). 

2. Observations and Technique 

2.1. The HESS Instrument 

HESS (the High Energy Stereoscopic System) is an array of 
four atmospheric Cherenkov telescopes located in the Khomas 
highland of Namibia at an altitude of 1800 m above sea-level. 
Each telescope consists of a 107m 2 optical reflector made up 
of segmented mirrors that focus light into a camera of 960 
photo-multiplier tube pixels (Bernlohr et al. 2003). The tele- 
scopes image the UV/blue flashes of Cherenkov light emitted 
by the secondary particles produced in gamma-ray-induced air- 
showers. Stereoscopic shower observations using the imaging 
atmospheric Cherenkov technique (e.g. Hillas 1996; Weekes 
1996; Daum et al. 1997) allow for accurate reconstruction of 
the direction and energy of the primary gamma rays as well as 
for the rejection of background events from air showers of cos- 
mic ray origin. HESS is sensitive to gamma rays above a post- 
cuts threshold energy of approximately 150 GeV and has an 
average energy resolution of ~ 16% (Aharonian et al. 2006b). 
Additionally, the high angular resolution (~0.1°), large field- 
of-view (~5°), and good off-axis sensitivity of the HESS array 
make it well suited for extended sources and scan-based obser- 
vations, where the source position is not known a priori. 

2.2. Data 

The observations discussed here were taken as part of the on- 
going HESS Galactic plane survey which currently covers the 
band -50° < / < 60° in galactic longitude and -3° < b < 3° 
in latitude. Data were taken as a series of 28-minute observa- 
tions (runs) centered on regular grid points covering the survey 
area. Additionally, several established sources were observed 
with pointed follow-up observations in wobble mode, where 
data are taken with an alternating offset from the target position 
of typically +0.7° in right ascension or declination. The set of 
usable runs were selected based on a standard set of hardware 
and weather conditions (Aharonian et al. 2006b). The sources 
in this study were chosen by selecting all locations in the HESS 
Galactic plane scan data set that have a pre-trials detection sig- 
nificance (with a fixed integration radius of 0.22°) greater than 
6cr (corresponding to a post-trials significance of 4cr, based on 
the very conservative estimate for the number of trials given in 
Aharonian et al. (2006d)), and for which no obvious cataloged 
counterpart can be associated (based on the criteria given in 
Section 2.4). Sources that were previously published (e.g. in 
Aharonian et al. 2006d) were excluded, except those that have 
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Source 


R.A 




Dec 


cr src (') 


HESS 11303-631$ 


13 h 03 m 


00 s 


-63°11'55" 


9.6 


HESS J1614-518$ 


16 h 14 m 


19 s 


-51°49'12" 


13.8 


HESS J1632-478 


16 h 32 m 


09 s 


-47°49'12" 


12.0 


HESS J1634-472 


16 h 34 m 


58 s 


-47° 16' 12" 


6.6 


HESS J1745-303 


17 h 45 m 


02 s 


30°22'12" 


12.6 


HESS J1837-069 


18 h 37 m 


38 s 


-6°57'00" 


7.2 


TeV 12032+4130$ 


20 h 32 m 


57 s 


41°29'57" 


6.2 



Table 2. Previously published unidentified VHE sources, not dis- 
cussed in this paper. Coordinates are in 12000 epoch, tr src is the 
intrinsic source extent (taking into account the instrumental re- 
sponse). Sources with $ have no obvious longer-wavelength counter- 
part. HESS 11632-478 has a possible HMXB counterpart, but the VHE 
source is extended; HESS 11634-472 may be related to an unidentified 
INTEGRAL source or nearby SNR, but is offset and morphologically 
dissimilar; HESS 11745-303 is partially coincident with an unidenti- 
fied EGRET source; and HESS 11837-069 is coincident with an as yet 
unidentified ASCA source. Results are from Aharonian et al. (2005c), 
Aharonian et al. (2006d), and Aharonian et al. (2005a). 



had increases in significance over 3 cr due to subsequent obser- 
vation. The eight sources that pass these selection criteria and 
their center positions (based on a model fit described in §2.3) 
are summarized in Table 1 . For reference, a summary of pub- 
lished results on previously reported unidentified VHE objects 
is given in Table 2. 

2.3. Analysis Technique 

The data presented here were analyzed using the standard 
HESS analysis scheme: calibrations are applied to the raw 
shower images (Aharonian et al. 2004) followed by an im- 
age cleaning procedure which removes noise due to fluctua- 
tions in the optical night-sky background light. The images are 
then parametrized using the Hillas moment-analysis technique 
(Hillas 1996), and gamma-ray selection criteria based on the 
image parameters are applied (Aharonian et al. 2006b). To re- 
duce systematic effects in the spectrum due to off-axis sensi- 
tivity that arise when images fall near the camera edge, an ad- 
ditional cut is applied to accept only data runs which are taken 
within an angular distance \j/ from the respective position of the 
object under analysis. For the spectral analysis, this is conser- 
vatively set to 2.0° to minimize systematic errors on the energy 
estimates (providing an average offset of 1.0° + 0.1°), while 
for the generation of the sky maps, it was set to 2.5° to maxi- 
mize the number of photons detected (giving an average offset 
of 1.9° ± 0.2°). Images from events passing the cuts for each 
telescope are combined to reconstruct the shower direction and 
energy. In the data presented here, two sets of gamma-ray selec- 
tion criteria are used to suppress events with hadronic origin: 
standard cuts, which are optimized using a simulated source 
with an energy spectrum with photon index T = 2.6 and a flux 
that is 10% of the Crab Nebula (a standard bright gamma-ray 
source) at VHE energies, and hard cuts which are optimized 
for a harder spectrum source (T = 2.0) with a flux that is 1% 
of the Crab Nebula. Standard cuts have an intrinsically lower 
energy threshold, but are looser and accept more background 



events, while the hard cuts provide better gamma-hadron sep- 
aration, and thus higher signal-to-noise ratio, at the expense of 
an increased energy threshold. Unless otherwise noted, hard 
cuts are employed for the spectral and morphological analy- 
ses presented in this article since they provide smaller system- 
atic errors due to a higher analysis energy threshold and better 
background rejection, though both sets are applied to check for 
consistency. 

The sky maps used for determining the source location and 
morphology are generated by accumulating the points of origin 
of each gamma-ray candidate in a two-dimensional histogram, 
subtracting a background map modeled by counting the num- 
ber of events which fall within an annulus (of average radius 
0.5°) about each grid point, excluding emission regions (the 
ring-background model described in Berge et al. 2007). The 
background is corrected for acceptance variations across the 
field of view. As an additional check, a background model us- 
ing the radial gamma-ray acceptance profile (as determined by 
dedicated off-source observations and simulations) in the field 
of view of each run is also used and compared for consistency. 
An elongated two-dimensional Gaussian convolved with the 
HESS point-spread function is fit to the resulting excess map 
to determine the centroid position, position angle, and extent 
of the source. To define the full extent of the source for spec- 
tral analysis, a histogram of the squared distance of each event 
to the fit position (6 2 ) is generated. The statistical significances 
of each excess measurement are calculated from the measured 
number of on- and off-source (background) events following 
the likelihood ratio procedure outlined in Li & Ma (1983). 

The background for spectra is estimated using the reflected- 
region technique where background events are selected from 
circular off-source regions in the field of view that have 
the same angular size and offset from the observation cen- 
ter position as the on-source region (Aharonian et al. 2006b). 
Background regions containing other known sources are ex- 
cluded. This technique provides a more accurate estimation of 
the background than the field-of-view model (described above) 
used to generate the sky maps, but is not as well suited for the 
generation of two-dimensional images. 

Spectra are generated following Aharonian et al. (2006b) 
for all events that fall within an angular distance 6{ al of the 
target position. This radius is chosen for each source as the 
distance where the radial excess distribution falls to a level in- 
distinguishable from noise (i.e. fully encloses the source). This 
provides a less biased estimate of the spectrum since it makes 
no assumption on the source morphology, but it decreases the 
signal-to-noise ratio since some additional background is in- 
cluded compared to an angular cut optimized for best signif- 
icance. An energy estimate for each event is calculated based 
on a comparison of the event's impact parameter, zenith angle, 
offset from the center of the field of view, and the amplitude 
of the integrated image for each telescope. The energy esti- 
mates for all events in the on and off-source regions are put 
into two histograms, which are then corrected for differing ex- 
posure, subtracted, and a flux is calculated for each energy bin 
by dividing by the observation time and the effective collection 
area of the telescopes (which is a function of energy, offset from 
the camera center, zenith angle, and the angle with respect to 
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Source 




Right Ascension 


Declination 


m 


bn 


Time (hrs) 


S (cr) 


Excess ( 


HESS J1427- 


-608 


14 h 27 ra 52 s 


-60°51'00" 


314.409 


-0.145 


21 


7.3 


197 


HESS J1626- 


-490 


16 h 26 m 04 s 


-49°05'13" 


334.772 


0.045 


12 


7.5 


153 


HESS J1702- 


-420f 


17 h 02 m 44 s 


-42°00'57" 


344.304 


-0.184 


9 


12.8 


412 


HESS J1708- 


-410t 


17 h 08 m 24 s 


-41°05'24" 


345.683 


-0.469 


39 


10.7 


513 


HESS J1731- 


-347 


17 h 31 ra 55 s 


-34°42'36" 


353.565 


-0.622 


14 


8.1 


218 


HESS J1841- 


-055 


18 h 40 m 55 s 


-05°33'00" 


26.795 


-0.197 


26 


10.6 


346 


HESS J1857- 


1-026 


18 h 57 m ll s 


02°40'00" 


35.972 


-0.056 


21 


8.7 


223 


HESS J1858- 


1-020 


18 h 58 m 20 s 


02°05'24" 


35.578 


-0.581 


25 


7.0 


168 



Table 1. Positions in equatorial (J2000 epoch) and Galactic (l,b) coordinates along with the detection significances of unidentified sources in 
the HESS Galactic Plane scan discussed in this paper. S is the significance (number of standard deviations above the background level) of the 
source using a fixed integration radius of 0.22°, which was used for selecting the sources from the scan data. The position of each source is 
based on a model fit to the background-subtracted gamma-ray maps (discussed in §2.3 and Table 3). The fit positions have an average statistical 
error of 0.05 degrees. Sources marked with a t are previously published in Aharonian et al. (2006d) and have been updated with new data. The 
exposure time is corrected for the off-axis sensitivity of the telescope system and accounts for instrumental readout dead-time. 



the Earth's geomagnetic field, as determined from simulations). 
The resulting fluxes are fit by a power-law of the form 

F(£) = ^i4vf (1) 

where F is the photon index and A^o is the flux normaliza- 
tion. Muon images are used to correct the energy estimate for 
changes in the optical efficiency of the telescopes over time 
(due to, e.g. the degradation of the mirrors) (Aharonian et al. 
2006b). The systematic error on the flux is conservatively esti- 
mated from simulated data to be 20% while the photon index 
has a typical systematic error of +0.2. 

To check the robustness of the results presented in this ar- 
ticle, the analysis has been repeated using several other back- 
ground models as well as with a completely separate analy- 
sis and calibration chain which used independent simulations 
and the forward-folded spectrum reconstruction technique de- 
scribed in Piron et al. (2001). 

2.4. Counterpart Search 

A search for counterparts to the VHE emission was made by 
first looking in source catalogs for objects which are of a type 
known to produce VHE photons, including the ATNF pulsar 
catalog (Manchester et al. 2005), the Green's supernova rem- 
nant catalog (Green 2004), and the High-Mass X-ray binary 
(HMXB) catalog by Liu et al. (2006). We also checked the 
Low-Mass X-ray binary (LMXB) catalog by Liu et al. (2007), 
the INTEGRAL source catalog (Bird et al. 2007), and the 
SIMBAD database. Sky maps for longer-wavelength survey 
data in the radio and X-ray wavebands, from the Molonglo 
(Green et al. 1999; Mauch et al. 2003), NRAO VLA (Condon 
et al. 1998), ROSAT (Voges et al. 2000), ASCA (Tanaka et al. 
1994) Galactic plane surveys, were compared with the HESS 
excess maps. Additionally, pointed observations made by the 
Chandra and XMM-Newton instruments were checked when 
available in the respective archives. Unless otherwise noted, 
ROSAT survey data between 1.0-2.4 keV and ASCA data be- 
tween 2-10 keV have been used. 

To reduce the number of chance coincidences with cata- 
loged sources, some loose selection criteria were applied: 



- Based on previous detections in the VHE energy range (e.g. 
Aharonian et al. 2006c,e), we consider the association of 
a VHE source with a shell-type SNR plausible only if the 
VHE emission roughly matches the angular size of the rem- 
nant and is not significantly offset. 

- Due to the large number of cataloged pulsars in the galac- 
tic plane, only those which are energetic enough to power 
a PWN which could produce VHE emission were consid- 
ered. A useful quantity for determining the possibility of 
VHE emission from pulsars is the spin-down flux measured 
at the solar-system, E/D 2 , where E is the spin-down lumi- 
nosity and D 2 is the distance to the object (both measur- 
able quantities) (Fierro et al. 1995). Defining the conver- 
sion efficiency, r\, as the ratio of the integral energy flux 
of a gamma-ray source over a typical energy range (e.g. 
200 GeV to 20 TeV) to the pulsar spin-down flux at the 
solar system, we find that for typical spectral character- 
istics of the sources discussed here, E/D 2 must be well 
above 10 33 ergsec -1 kpc 2 to produce the observed emis- 
sion, even assuming 100% efficiency; for this reason, pul- 
sars with lower spin-down fluxes are not plotted in the fig- 
ures given later in this paper. In cases where the distance 
estimate is not known, we assume a distance of 3 kpc. We 
note that efficiencies greater than 100% are not completely 
excluded, since the spin-down flux might have been higher 
in the past and the particle cooling times might be compa- 
rable to the pulsar's age (Aharonian et al. 2007d). However, 
to claim a plausible identification of a VHE source with a 
pulsar, we require efficiencies of < 10% and a reasonably 
small angular distance for the purpose of this study to keep 
the number of chance coincidences low, unless there are 
other multi-frequency data that would support the associa- 
tion. 

- XRBs from which VHE emission is established are 
HMXBs that either exhibit a jet (e.g. LS 5039, Aharonian 
et al. 2006f), or where the compact object is a pulsar pow- 
ering a PWN (e.g. PSR B1259, Aharonian et al. 2005d); all 
appear variable and point-like in the VHE band. We believe 
the chance probability of the appearance of an XRB within 
the 3cr contours of an extended HESS source to be reason- 
ably low, therefore we discuss such associations, but for the 
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Source 


<T\ ( ) 


CT2 ( ) 


Angle ( ) 


HESS J1427-608 


0.04 ± 0.02 


0.08 ± 0.03 


80 ± 17 


HESS J1626-490 


0.07 ± 0.02 


0.10 ±0.05 


3 ± 40 


HESS J1702-420 


0.30 ± 0.02 


0.15 ±0.01 


68 ± 7 


HESS J1708-410 


0.06 ± 0.01 


0.08 ± 0.01 


-20 ± 23 


HESS J1731-347 


0.18 + 0.07 


0.11 ±0.03 


-89 ± 21 


HESS J1841-055 


0.41+0.04 


0.25 ± 0.02 


39 ± 6 


HESS J1857+026 


0.11 ±0.08 


0.08 ± 0.03 


-3 ± 49 


HESS J1858+020 


0.08 ± 0.02 


0.02 ± 0.04 


4 ± 17 



Table 3. Results from an elongated 2-D Gaussian model fit (see §2.3) 
to the gamma-ray excess for each source. cr l and <x 2 are the intrinsic 
semi-major and semi-minor axes (in degrees on the sky), with the ef- 
fect of the point-spread function removed. The errors are statistical. 
The position angle is measured counter-clockwise in degrees relative 
to the RA axis. 



moment ignore XRBs lying outside the sources. Although 
the possibility exists that such an object might also power 
an extended, possibly asymmetric VHE source (e.g. Cheng 
et al. 2006), this has so far not been observed. 

3. Results 

Results of the size and spectral fits for each source are sum- 
marized in Tables 3 and 4, respectively. The spectrum for each 
source is plotted in Figure 8. In the following sections, a de- 
tailed discussion of each source and related cataloged sources 
or hot-spots within each field of view is given. 

3.1. HESS J1 427-608 

HESS J1427-608 (Figure 1) is located approximately 1° away 
from the hard X-ray and GeV gamma-ray source G3 13.2+0.3 
(a strong radio source located in the Kookaburra complex) 
(Aharonian et al. 2006a), and has a slightly extended morphol- 
ogy consistent with a symmetric Gaussian of radius cr = 3'. Its 
spectrum is fit by a power-law with index 2.2 + 0.1 stat + 0.2 sys . 
Radio and X-ray survey data of the region (overlaid in Figure 
1 from the Molonglo and ROSAT surveys, respectively) show 
no evidence for significant emission at distances of 0.5° or 
closer to the centroid position of HESS J1427-608. There are 
no nearby pulsars or supernova remnants, and an association 
of HESS J1427-608 with the unidentified INTEGRAL source 
IGR J14331-61 12 is unlikely due to the large angular distance 
separating the two sources. 

3.2. HESS J1 626-490 

HESS J1626-490, located exactly on the Galactic plane 
(Figure 2), is a gamma-ray source with an approximately 
radially-symmetric Gaussian morphology (with 5' extent), and 
a power-law energy spectrum with photon index 2.2 + 0.1 stat + 
0.2 sys . There is a slight extension toward increasing right ascen- 
sion which is only marginally significant, but may be an indica- 
tion of a second VHE source. Within the gamma-ray emission 
region, there exists some weak radio emission, along with the 
unidentified X-ray source 1RXS J162504-490918, which lies 



approximately 10' from the centroid position and is a possible 
X-ray counterpart. This X-ray source, marked with an "X" in 
the figure, has an extent of 13" and an absorption-corrected flux 
between 0.1-2.0 keV of 1.7 x 10~ 13 erg ctrT 2 s _1 , assuming a 
photon index of 2.0 (Voges et al. 2000; Mukai 1993). The shell- 
type supernova remnant G335.2+00.1 (MSH 16-44) (Whiteoak 
& Green 1996) lies just outside the significant emission region 
of HESS J1626-490, as does the LMXB 4U 1624-490 (Smale 
et al. 2000), and the HMXB IGR 16283-4838 (Bird et al. 2007), 
which are not considered plausible candidates due to their off- 
sets. 

3.3. HESS J1 702-420 

First discovered by HESS at an approximately 6<x significance 
level (Aharonian et al. 2006d), HESS J1702-420 (Figure 3) 
is now seen with increased observation time at a significance 
level of 13<x. Its spectrum is characterized by a power-law with 
index 2.1 + 0. l stat ± 0.2 sys , slightly harder than the previously 
reported value of 2.3 +0.15 sta t±0.2 sys , which was derived from 
a smaller integration radius, less statistics, and over a smaller 
energy range. The results, including the source location, are 
consistent within the errors. The emission "tail" extending to 
positive galactic longitude and latitude is statistically signifi- 
cant, giving the source an elongated morphology (see Table 3). 
The nearby pulsar PSR J 1702-4 128 (to the north of the VHE 
emission region, Figure 3) lies at the edge of the gamma-ray 
emission, and with E/D 2 = 1.3 • 10 34 erg s _1 kpc~ 2 , it provides 
enough spin-down energy loss to produce the observed emis- 
sion (assuming a rather high conversion efficiency of ~ 70% if 
the present distance estimate of 5 kpc is correct) and may be a 
counterpart if it powers an extremely asymmetric pulsar wind 
nebula. The nearby shell-type supernova remnant G344. 7-00.1 
(seen in the radio image) is also detected by ASCA in the 2-10 
keV X-ray energy band (Sugizaki et al. 2001), however is an 
unlikely counterpart due to its small angular size and distance 
from the peak of the emission region. Three X-ray binaries are 
also located nearby the source (see the figure), but are outside 
the significant emission region. 

3.4. HESSJ1708-410 

HESS J1708-410 (Figure 4), situated between the super- 
nova remnant RXJ 1713.7-3946 (Aharonian et al. 2006e) and 
HESS J1702-420, was first reported at a significance level of 
approximately 7<x (Aharonian et al. 2006d). With additional 
observations of the region (mostly from the edge of pointed ob- 
servations centered on RXJ1713. 7-3946), the data set now has 
a statistical significance of 1 l<x. The spectrum is fit by a power- 
law with index 2.5 ± 0. 1 stat + 0.2 sys , which is slightly softer than 
the previously published result of 2.3 ± 0. 1 stat ± 0.2 sys made with 
lower statistics, a smaller integration radius, and over a smaller 
energy range (Aharonian et al. 2006d), though is within er- 
rors. The compact morphology of HESS J1708-410 is consis- 
tent with a slightly elongated Gaussian of approximately 0.08° 
extent, with no significant emission beyond 0.3°, ruling out 
SNR G345.7-00.2 or nearby radio hot-spots as obvious coun- 



6 



The HESS Collaboration: Unidentified HESS Sources 



Source 


&int 


Time 


S 


Excess 


r 


N ■ 10~ 12 




^Vnax 


D.O.F. 




V ) 


(first 


(rrt 


( POllTlf 5^ 




Ci-in - 2 C — 1 Tf=V — 1 t 

(cm s lev ) 




(TeV 1 ) 




HESS J1427-608 


0.2 


21 


5.6 


165 


2.16 + 0.14 


1.3 ±0.4 


0.97 


50 


4.6/6 


HESS J1626-490 


0.5 


8 


6.0 


167 


2.18 ±0.12 


4.9 ± 0.9 


0.60 


50 


6.4/7 


HESS J1702-420 


0.6 


7 


10.6 


596 


2.07 ± 0.08 


9.1 ± 1.1 


0.50 


50 


11.3/7 


HESS J1708-410 


0.3 


45 


10.3 


542 


2.46 ± 0.08 


2.7 ± 0.3 


0.50 


60 


5.4/5 


HESS J1731-347 


0.6 


11 


8.3 


495 


2.26 ±0.10 


6.1 ±0.8 


0.50 


80 


2.8/6 


HESS J1841-055 


0.7 


10 


10.7 


723 


2.41 ±0.08 


12.8 ± 1.3 


0.54 


80 


10.4/6 


HESS J1857+026 


0.46 


15 


10.2 


425 


2.39 ± 0.08 


6.1 ±0.7 


0.60 


80 


5.9/6 


HESS J1858+020 


0.15 


23 


7.4 


116 


2.17 ±0.12 


0.6 ±0.1 


0.50 


80 


4.9/6 



Table 4. Summary of spectral parameters for each source from a power-law fit to the spectral data (E = NoE~ T ) over the energy range 
£min _ £ ma i- The integration radius, 6 inl is chosen to fully enclose each source. Only data with observation positions offset less than 2° from the 
source position were included. The errors shown are statistical; the systematic error is conservatively estimated to be 20% on the flux and ±0.2 
on the spectral index. Plots of all spectra are given in Figure 8. 



terpart candidates. Although several ROSAT hard-band X-ray 
hot spots exist in the field-of-view (e.g. the XRB 4U 1708-40 or 
1RXS J17101 1.5-405356, see figure), the closest is 0.2° away 
and is not obviously connected with the gamma-ray emission. 
There is an XMM-Newton exposure centered on G345.7-00.2, 
in which no significant emission is seen near the VHE position. 
Additionally, an ASCA exposure of the region reveals only a 
single point-like source located over a degree from the HESS 
source. 

3.5. HESSJ1731-347 

HESS J1731-347 (Figure 5) is detected at an ~8 cr level, ex- 
hibiting a power-law spectral index of 2.3 + 0.1 sta t ± 0.2 sys . The 
source has a significant tail which extends westward, giving it a 
non-Gaussian morphology, possibly indicating the presence of 
more than one or an extended non-uniform source. A slice in 
the uncorrelated excess event map along the axis of the emis- 
sion does not show a conclusive separation between the two 
"peaks", and a spectral analysis of each gives the same photon 
index within systematic errors. For this reason, the emission is 
treated here as a single source. 

A bright X-ray point source (1RXS J173030.3-343219, 
labeled as "X" in the figure) is seen in the ROSAT data, 
approximately 0.4 degrees in the direction of the Galactic 
Plane from the centroid position, and has an absorption- 
corrected flux in the 0.1-2.4 keV range of approximately 
2.0 x 10" 11 erg cirT 2 sec" 1 (Voges et al. 1999; Mukai 1993), 
assuming a spectral index of 2.0. This source is identified 
with the cataclysmic variable (CV) star HD 158394, and is 
not expected to produce VHE emission. However, around the 
brightest part of the TeV emission, there is some unidenti- 
fied nebular X-ray emission that partially matches the mor- 
phology of the HESS source, and may well be the X-ray 
counterpart. This diffuse X-ray emission includes the extended 
ROSAT source 1RXS J 17325 1.1 -344728 (labeled X x in the 
figure), which has an extension of 2' and X-ray flux of (7 + 
1) x 10~ 12 erg crrT 2 sec -1 , and a nearly coincident point-like 
radio source labeled 353.464-0.69 in the VLA survey data 
(Zoonematkermani et al. 1990); their association with the VHE 
emission is unclear. The strong point-like radio source 173028- 
344144 (Condon et al. 1998), labeled R in the figure, also 



lies to the right of the peak of the VHE emission. The X-ray 
emission about a degree away to the north in the figure comes 
from the LMXB GX 354-0, however due to its distance from 
HESS J1731-347 and since these objects are not known to 
produce extended gamma-ray emission, it is an unlikely coun- 
terpart candidate. No known high spin-down flux pulsars lie 
within the emission region. 

3.6. HESSJ1841-055 

HESS J1841-055 (Figure 6) exhibits a highly extended, possi- 
bly two or three-peaked , morphology; however, the "dip" be- 
tween the peaks along the major axis is not statistically signifi- 
cant (< 1.5cr). The source has a spectrum that is fit by a power 
law with index 2.4 + 0.1 stat + 0.2 sys . An association with ei- 
ther pulsar PSR J1841-0524 (E/D 2 = 4.4 • 10 33 erg s 1 kpc 2 ) 
or PSR J1838-0549 (E/D 2 = 4.7 • 10 33 erg s 1 kpc 2 ), is not 
ruled out, however taken separately, each would require ap- 
proximately 200% efficiency to explain the VHE emission. 
This is not completely implausible if both pulsars contribute 
together or if either had a much higher spin-down luminosity 
in the past (PSR J 1838-0549 is estimated to have a relatively 
old characteristic age of 112 kyr, while PSR J 184 1-0524 is 
about 30 kyr old (Manchester et al. 2005)). PSR J1837-0604 
(E/D 2 — 5.2 • 10 34 erg s _1 kpc~ 2 ) has a high enough spin-down 
flux to be a counterpart candidate, however it is well outside the 
emission region. There are no cataloged PWN at longer wave- 
lengths identified with any of the three pulsars (e.g. Gotthelf 
2004). The SNR G027.4+00.0 (also known as Kes 73), which 
is visible in both X-ray and radio wave bands, lies at the edge 
of the emission, though does not appear related due to its small 
angular size. Additionally, the high-mass X-ray binary J 1839- 
06 also lies near the edge of the significant TeV excess. 

ASCA observations of the Scutum arm region reveal a 
point-like source, AX J1841. 0-0536, near the center of the 
VHE emission, which based on its X-ray light curve and op- 
tical emission is suggested to be a Be/X-ray binary pulsar 
(Bamba et al. 2001) with a flux in the 6-20 keV energy range of 
l.lxlO~ 10 ergcrrT 2 s _1 and photon index of 2.2+0.3 (Filippova 
et al. 2005). A Chandra observation of this object confirms the 
identification, with a flux in the 0.5-10 keV energy range of 
4.2 x 10~ 12 erg cm" 2 s" 1 (Halpern et al. 2004). Given its point- 
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like extent, AX J1841.0-0536 is not large enough to explain the 
entire HESS source, however it may well be responsible for a 
component of the emission. 

Also within the VHE emission region lies the diffuse source 
G26.6-0.1, which was detected in the ASCA Galactic Plane 
Survey and is postulated based on its spectrum to be a candi- 
date supernova remnant (Bamba et al. 2003), and is also coin- 
cident with an H II region (Lockman 1989). With its angular 
size of 8.3' (FWHM), small distance (approximately 1.3 kpc), 
and non-thermal spectrum, this object also may also contribute 
to a component of the VHE emission. Additionally, the nearby 
source AX J18406-0539 is possibly an XRB at a distance of 
1.1 kpc (Masetti et al. 2006), though given positional errors, 
may well be the same source as AX J1841.0-0536 (Negueruela 
& Schurch 2007). 



3.7. HESS J1 857+026 

HESS J 1857+026 (Figure 7) is an approximately radially- 
symmetric extended VHE gamma-ray source located on the 
Galactic Plane. The source is detected by HESS at a 9<x sig- 
nificance level at energies above 300 GeV and has a differen- 
tial spectral index of 2.4 + 0.1 stat ± 0.2 sys . The slight extension 
of the source seen toward the north is significant (~ 5<x) and 
may indicate a more extended morphology or the presence of 
a weaker nearby source, though more observation time will be 
needed to make a conclusive statement. 

This source lies approximately 0.7° from HESS J1858+020 
(see §3.8), which is most probably a separate source since no 
significant emission connects the two. An association with the 
supernova remnant G036.6-00.7, which lies over a degree from 
the centroid position, is unlikely. Though an ASCA observa- 
tion exists which is roughly centered on the source position, 
no excess was seen, implying a 95% absorbed flux upper- 
limit of 1.2 • 10~ 12 ergcnr 2 s" 1 between 2-10 keV The X- 
ray source seen about a quarter of a degree from the cen- 
troid position is the point-source 1 RXS J 1 85609.2+02 1 744 (la- 
beled X in the Figure, and coincident with the ASCA source 
AXJ 185608+0218), which has a flux in the 0.1-2.4 keV range 
of (0.32 ± 0.06) x 10~ 12 erg crrT 2 s _1 , assuming a photon in- 
dex of 2.0; its distance from the emission region makes it an 
unlikely counterpart candidate, however. 



3.8. HESS J1 858+020 

The weak gamma-ray source HESS J1858+020 (shown also in 
Figure 7) lies close to HESS J 1857+026; however, there is no 
significant emission connecting them, suggesting that they are 
distinct objects. It is detected at a significance level of 7<x with 
a differential spectral index of 2.2 ± 0.1 stet ± 0.2 sys . Though 
nearly point-like, its morphology shows a slight extension of 
~ 5' along its major axis. PSR J1857+0143 (E/D 2 = 1.7 • 
10 34 ergs _1 kpc -2 ) is powerful enough to explain the source, 
but is significantly offset. 



4. Summary 

The eight VHE gamma-ray sources discussed here are all ex- 
tended objects with angular sizes ranging from approximately 3 
to 18 arc minutes, lying close to the Galactic plane (suggesting 
they are located within the Galaxy). In each case, the spectrum 
of the sources in the TeV energy range can be characterized as 
a power-law with a differential spectral index in the range 2.1 
to 2.5. The general characteristics of these sources — spectra, 
size, and position — are similar to previously identified galac- 
tic VHE sources (e.g. PWNe), however since these sources 
have so far no clear counterpart in lower-energy wavebands, 
further multi-wavelength study is required to understand the 
emission mechanisms powering them, and therefore follow-up 
observations with higher-sensitivity X-ray and GeV gamma- 
ray telescopes will be beneficial. Since most VHE sources are 
predicted to emit X-ray and radio emission, a non-detection 
of longer-wavelength emission with current-generation exper- 
iments for some of these objects may be an indication that a 
new VHE source class exists (as suggested by Aharonian et al. 
2005b), and may provide new insight into high-energy pro- 
cesses within our Galaxy. 
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Fig. 1. Left: A VHE gamma-ray image of HESS J1427-608 (center position marked with a star). The image is of gamma-ray excess counts 
smoothed with a Gaussian filter with standard deviation 0.1° (referred to as the smoothing radius hereafter). The smoothing radius is chosen 
according to event statistics and therefore differs from source to source and is shown as a green circle in the lower left-hand corner. The color 
scale of the image is set such that the blue/red transition occurs at approximately the 3cr (pre-trials) significance level. Overlaid on the image are 
the significance contours starting at 4cr in lcr steps. The Galactic plane is marked with a dashed line. The gamma-ray excess at the right of the 
image is the hard X-ray source known as the Kookaburra/Rabbit, which is discussed in Aharonian et al. (2006a). Right: The HESS significance 
contours (black) overlaid on a radio image (Green et al. 1999) (grey-scale, in Jy/beam). The green contours are from a ROSAT hard-band X-ray 
image (Voges et al. 2000) which has been adaptively smoothed with the FTOOLS /adapt algorithm to accentuate diffuse emission (Blackburn 
1995). Also plotted are ATNF pulsars with E/D 2 > 10 33 erg S" 1 kpc 2 , SNRs from Green's catalog, HMXBs and LMXBs from the catalogs of 
Liu et al, and INTEGRAL sources (see Section 2.4 for references). In this case, only the INTEGRAL source IGR J14331-61 12 lies within the 
field of view. 
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Fig. 2. Left: A VHE gamma-ray image of HESS J 1626-490 plotted as in Figure 1, with a smoothing radius of 0.1°. Right: the HESS significance 
(black) and adaptively smoothed ROSAT X-ray contours (green), overlaid on the Molonglo radio image (grey-scale). Also plotted is the SNR 
G335.2+00.1 (circle marking extent), the HMXB IGR 16283-4838, the LMXB 4U 1624-490, and the unidentified X-ray source 1RXS J162504- 
490918 (labeled X). 
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Fig. 3. Left: A VHE gamma-ray image of HESS J1702-420, plotted as in Figure 1, with a smoothing radius of 0.06°. Right: the HESS 
significance (black) and adaptively smoothed ROSAT X-ray contours (green), overlaid on the Molonglo radio image (grey-scale). Also plotted 
are the positions of the SNR G344.7-00.1 (circle), three HMXBs (magenta squares), and the high spin-down flux pulsar PSR J1702-4128 (red 
triangle). 
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Fig. 4. Left: A VHE gamma-ray image of HESS J1708-410 plotted as in Figure 1, with a smoothing radius of 0.06°. The slight excess seen 
on the lower-right corner of the image is HESS J1702-420 (see previous section), while the excess seen at the upper-left corner is part of 
RXJ 1713.7-3946 (Aharonian et al. 2006e). Right: The HESS significance (black) and adaptively smoothed ROSAT X-ray contours (green), 
overlaid on the Molonglo radio image (grey-scale). Also plotted are the positions of the SNR G345 .7-00.2, the ROSAT source 1RXS J171011.5- 
405356 (labeled "X"), and the LMXB 4U 1708-40 (square). 
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Fig. 5. Left: A VHE gamma-ray image of HESS J1731— 347 plotted as in Figure 1, with a smoothing radius of 0.1°. Right: The HESS signifi- 
cance (black) and adaptively smoothed ROSAT X-ray contours (green), overlaid on the Molonglo radio image (grey-scale). Also shown are the 
positions of a high spin-down flux pulsar (filled triangle), the low-mass X-ray binary GX 354-0, and cataloged supernova remnants (blue circles 
marking extent). The source labeled Xi is the ROSAT source 1RXS J17325 1.1-344728, R 2 is the point-like radio source 173028-344144, and 
X 2 is the ROSAT point-source 1RXS J173030.3-343219. 
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Fig. 6. Left: A VHE gamma-ray image of HESS J1841-055, plotted as in Figure 1, with a smoothing radius of 0.07°. Right: The HESS 
significance (black) and adaptively smoothed ROSAT X-ray contours (green), overlaid on the NVSS radio image (grey-scale). Also shown are 
the positions of known high spin-down flux pulsars (filled triangles), the SNR Kes 73 (circle), the X-ray candidate SNR G26.6-0. 1 (blue circle), 
and the HMXB AXJ 1841.0-0536. The Ginga source GS 1839-06 is compatible with the location of AXJ 1841.0-0536. 
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Fig. 7. Left: A VHE gamma-ray image of HESS J1857+026 and HESS J1858+020, plotted as in Figure 1, with a smoothing radius of 0.08°. 
Right: The HESS significance (black) and adaptively smoothed ROSAT X-ray contours (green), overlaid on an NVSS radio image (grey-scale). 
Also shown are the positions of a known high spin-down flux pulsar (filled triangle), the SNR G036. 6-00.7 (circle), the ROSAT point-source 
1RXS J185609.2+021744 (labeled X), and the HMXB XTE J1858+034 (square). 
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Fig. 8. Spectra for each unidentified source, with power-law fits. See Table 4 for detailed fit information. 



